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(Dated: February 3, 2008) 

We present a search for a neutral particle, pair-produced in pp collisions at v / s=1.96 TeV, which 
decays into two muons and lives long enough to travel at least 5 cm before decaying. The analysis 
uses ^380 pb _1 of data recorded with the DO detector. The background is estimated to be about 
one event. No candidates are observed, and limits are set on the pair production cross section times 
branching fraction into dimuons + X for such particles. For a mass of 10 GeV and lifetime of 
4 x 10~ n s, we exclude values greater than 0.14 pb (95% C.L.). These results are used to limit the 
interpretation of NuTeV's excess of di-muon events. 



PACS numbers: 14.80 Ly, 12.60 Jv, 13.85 Rm 

Several models including supersymmetry with im- 
parity violation 0, Q and hidden valley theories Q pre- 
dict the existence of neutral, long-lived particles that give 
rise to a distinctive signature of two leptons arising from 
a highly displaced vertex. The Fermilab neutrino exper- 
iment NuTeV observed an excess of di-muon events that 
could be interpreted as such a signal 0, 0, 01 ■ Experi- 
ments at the CERN e+e" Collider (LEP) have looked for 
short-lived neutralino and chargino decays and longer- 
lived charged particles but did not search for this 
topology. 

In this Letter we present a search for a light, neutral, 
long-lived particle (N® L ) pair-produced in pp collisions at 
= 1-96 TeV and recorded with the DO detector, using 
380 pb _1 of data from Run II of the Fermilab Tevatron 
Collider. The final state under study is the decay of an 
N^l into two muons and possibly a neutrino after the 
Nj jL has traveled at least 5 cm. The particle is assumed 
to have a mass as low as several GeV. The analysis re- 
ported here explores a region of phase space previously 
unexplored by collider experiments. 

We use i?-parity violating (RPV) decays of neutrali- 
nos (x?) to [i + [i~v (Fig. [TJ as a benchmark model and 
to determine signal efficiency. Here the RPV couplings 
are expected to be small and lead to long lifetimes 
The results are applicable to any pair-produced neutral 
particle with similar kinematics. 

The DO detector consists of a central-tracking system, 
a liquid-argon and uranium calorimeter, and an outer 
muon system |l0|. Each of these is used in this analy- 
sis, with an emphasis on the muon system for particle 
identification and on the tracking system for momentum 
measurement and vertexing. 

The central tracker consists of a silicon microstrip 
tracker (SMT) and a central fiber tracker (CFT), both 
located within a 2 T superconducting solenoidal magnet. 
It is optimized for tracking and vertexing at pseudora- 
pidities \r\\ < 3 and \v\ < 2.5, respectively, where r\ = 
— ln[tan(0/2)] and 9 is the polar angle with respect to 
the proton beam direction. The CFT has 8 axial and 
8 stereo layers with an innermost (outermost) radius of 
20 (52) cm. The calorimeter consists of a central section 
(CC) covering \rj\ < 1.1, and two end calorimeters (EC) 
that extend coverage to M ss 4.2, with all three housed 
in separate cryostats |ll| . The outer muon system, at 
1 77 1 < 2, consists of a layer of tracking detectors and scin- 
tillation trigger counters in front of 1.8 T iron toroids, 
followed by two similar layers after the toroids [l2j| . 

We use the volume inside the CFT inner radius as a 
decay region. This allows the full CFT and muon sys- 
tems to be used for detection of decay products, ensur- 
ing robust track reconstruction and muon identification. 
Events are required to pass a di-muon trigger. 

The strategy is to identify events with at least two 
opposite-sign, isolated muons, defined as hits in the muon 




FIG. 1: Feynman diagrams for pair production (left) and de- 
cay (right) of a neutral particle, in this case neutralinos with 
_R-parity violation. 

system matched to a track in the CFT. Each pair is fit 
to a vertex. The signal sample uses events with muon 
vertices that are displaced more than 5 cm (in the plane 
transverse to the beamline) from the primary vertex. To 
characterize the displacement, we define the vertex radius 



r = y/(X- X PV f + (Y- Y PV f (1) 

where X ,Y are the x,y positions of the fit di-muon vertex 
and Xpv,Ypv are the x,y positions of the primary vertex 
(PV). DO uses a right-handed coordinate system with 
the positive z-axis defined by the proton direction and 
positive y-axis pointed upward. 

Studies of K s mesons are performed to test the re- 
construction efficiency for highly displaced vertices. We 
search for K s mesons in data and Monte Carlo (MC) 
by fitting track pairs to a common vertex and selecting 
those with an invariant mass around the K s peak. We 
are able to observe decay lengths greater than 20 cm 
and demonstrate that the data and MC follow the same 
radial dependence. The efficiency varies by 30% in the 
range r=5-20 cm. 

Selection criteria are chosen to minimize background 
while maintaining signal efficiency. All possible primary 
and secondary vertices are determined using tracks, ex- 
cept those associated with muons. The hard scatter ver- 
tex is determined by clustering tracks into seed vertices 
by a Kalman filter algorithm 13]. A probability function 
based on the pr of tracks attached to each vertex is used 
to rank the likelihood that it comes from a minimum 
bias interaction. The primary vertex is the one with the 
lowest probability. The PV is required to be within 0.3 
cm of the beamline in x and y and within 60 cm of the 
detector center in z. 

We require two muons which have hits in each of the 
three layers of the muon system, are matched to a track 
in the central tracker, have a good track fit, at least 14 
CFT hits associated with the track, transverse momen- 
tum > 10 GeV, and are isolated. Two methods are used 
to define isolation. First, the direction of the muon is 
projected to the calorimeter and the transverse energy in 
all cells within an annular cone 0.1<i?<0.4is summed 
(calorimeter isolation), where R = <J (A0) 2 + (A77) 2 . 
Second, the transverse momentum of all tracks within a 
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cone of R < 0.5 (except for the muon tracks) is summed 
(track isolation). Both the calorimeter and track isola- 
tions are required to be less than 2.5 GeV. Cosmic ray 
muons are rejected by requiring the time measured by 
the muon scintillator counters to be that expected for a 
particle produced at the nominal beam crossing time. To 
enhance the signal, both muons must have a distance of 
closest approach (DCA) of greater than 0.01 cm in the 
x — y plane and more than 0.1 cm along the z-axis from 
any vertices. The two muons must have an opening angle 
less than 0.5 radians and have opposite charge. 

All pairs of muons passing the above quality cuts are fit 
to a common vertex requiring a x 2 ' /N^ f < 4. The radial 
distance between the di-muon vertex and the primary 
vertex must be six times the resolution of the di-muon 
vertex measurement and be between 5 and 20 cm. This 
defines our signal region. 

We use data to estimate the background for this search. 
By allowing events to pass or fail two different selection 
criteria (the DCA and vertex radius cuts) we define four 
regions. For the DCA cut, we require either: (1) one 
track to pass the DCA cut and one to fail it, or (2) both 
tracks to pass the DCA cut. For the vertex radius we 
define two regions: (A) 0.3 < r < 5 cm, or (B) 5 < r < 
20 cm. This defines Samples 1A, 2A, IB, and 2B. 

We observe four events in Sample 1A, one event in 
Sample IB, and three events in Sample 2A. Assuming 
no correlation between selection criteria and no signal, 
the ratio of the number of events in region to 2B to the 
number in IB should equal the ratio of the number of 
events in region 2A to the number in 1A. This can be re- 
expressed to give an estimate of the background in the 
signal sample (Sample 2B): 



An 



bkgd 



Sample 2A 
Sample 1A 



x Sample IB. = 0.75 ± 1.1 events 

Due to bias from a correlation between the vertex ra- 
dius and DCA criteria, we perform tests of this method 
and assign a systematic uncertainty using several addi- 
tional samples. The spread in these results is used to 
assign a systematic uncertainty to account for the corre- 
lation between the vertex radius and DCA cut. Thus, we 
estimate the background in the signal region to be 0.75 
±1.1 (stat) ±1.1 (syst) events. 

Figurc|21shows the vertex radius distribution for events 
where one or both muons pass the DCA criteria. Exam- 
ination of the signal region yields events passing all 
criteria. Therefore we set a limit on the cross section as 
a function of lifetime. The lifetime dependence is cal- 
culated based on the fraction of events, /, which decay 
within our signal region. 

Signal MC events are generated using SUSYGEN 
and an unconstrained minimal supersymmetric model 
with R-parity violation Q, 0, pi using CTEQ5L parton 
distribution functions (PDFs) [l5|. The following param- 
eters are used: tan/3 = 10, = -5000, M 2 = 200 GeV, 
M 3 = 400 GeV, M squark = 300 GeV, M slcpton = M snu = 
Msbottom = Mstop = 1500 GeV. The Xi mass is about 
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FIG. 2: Distribution of the vertex radius for events where one 
muon passes the DCA criteria and the second fails it (top) and 
where both muons pass the DCA criteria (bottom). 



equal to the Mi parameter. Similar sets are generated 
with Mi = 3, 5, 8, 10, 15, 20, 30, and 40 GeV yielding 
pair production cross sections in the range 0.025-0.013 
pb. While the parameters are different than those used 
in other neutralino searches , they are chosen to give a 
model that provides a particle consistent with the NuTeV 
result and does not violate limits from LEP searches (in- 
cluding the Z boson invisible width). The lifetime is 
determined primarily by the slepton mass and the A122 
parameter, however, in the detector simulation we choose 
to ignore the lifetime and force exactly one of the two XiS 
to decay within a cylinder of radius 25 cm. The vertex 
is selected along the Xi trajectory such that the radius 
distribution is flat over the range 0-25 cm. The other Xi 
is required to escape the detector. The dependence of the 
acceptance on the lifetime is accounted for in the inter- 
pretation of the final result. The average x? transverse 
momentum (pt) is ~85 GeV. 

Our uncertainty estimate on the luminosity times sig- 
nal acceptance is summarized in Table H] The MC accep- 
tance uncertainty is statistical. Tracking, isolation, and 
muon reconstruction data/MC corrections are estimated 
using the Z boson mass peak, yielding 0.72 ± 0.07. The 
vertex reconstruction data/MC is found using K s events 
(0.92 ± 0.14). A PDF uncertainty on the signal efficiency 
of ±4% is assigned using the CTEQ6.1M PDF set 16]. 

These event numbers, efficiencies, acceptances and un- 
certainties are combined to set a 95% (99%) confidence 
level limit on the cross section a(pp — > Nl L N® L X) times 
branching fraction BF(A° L — > + X) as a function 

of the lifetime (Fig. using a Bayesian technique 0J 
and assuming zero background. The limit for a 10 GeV 

11 s is 0.14 pb (95% C.L.). 
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TABLE I: Acceptance, error, and limits for the MC sig- 
nal points. The lifetime acceptance is the factor by which 
the limit is adjusted due to the fraction of events which de- 
cay within the 5-20 cm region and is given for a lifetime of 
4 x 10 -11 s. The luminosity x acceptance includes the MC 
signal acceptance, the trigger efficiency, the data/MC correc- 
tion factors, and the luminosity. The limits are given for the 
same lifetime. 









Luminosity 


95% C.L. 


M(xl) 


Monte Carlo 


Lifetime 


x acceptance 


limit 


(GeV) 


acceptance 


acceptance 


(Pb- 1 ) 


(Pb) 


3 


0.095 ± 0.005 


0.51 


23.9 ± 4.8 


0.28 


5 


0.114 ± 0.005 


0.61 


28.7 ± 5.8 


0.19 


8 


0.141 ± 0.006 


0.67 


35.5 ± 7.1 


0.14 


10 


0.136 ± 0.006 


0.68 


34.3 ± 6.8 


0.14 


15 


0.139 ± 0.006 


0.65 


35.1 ± 7.0 


0.15 


20 


0.130 ± 0.005 


0.62 


32.8 ± 6.5 


0.17 


30 


0.099 ± 0.003 


0.55 


24.8 ± 4.9 


0.25 


40 


0.079 ± 0.004 


0.48 


20.0 ± 4.1 


0.35 




L = 380±25 pb 

D0 95% Exclusion 

| | D0 99% Exclusion 



i<r icr 1 1 
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FIG. 3: Limit on cross section x branching fraction for the 
pair-production of neutral, long-lived particles as a function 
of lifetime. The dark gray area and above represents the DO 
99% C.L. limit for the 5 GeV mass point. The solid line shows 
the DO 95% C.L. limit. The light gray region represents the 
NuTeV 99% C.L. exclusion |4| converted to a pp cross section 
at ,/s — 1960 GeV. The white region represents a 99% C.L. 
preferred region given the three events from NuTeV. 



Figure 0] shows how the DO limit varies with mass at a 
lifetime of 4 x lCT 11 s. 

In order to compare with DO, we convert the NuTeV re- 
sult from pp production at -^/s = 38 GeV to pp production 
at y/s — 1960 GeV using the ratio of cross sections for 
SUSY neutralino pair production calculated with the pa- 
rameters from our 5 GeV signal simulation. The NuTeV 
lifetime is converted from kilometers to seconds assum- 



ing an average momentum (along the neutrino beam di- 
rection) of 121 GeV. Given that NuTeV observed three 
events jj], a preferred region is found using the ratio of 
the 99% CL lower and upper limits on three events de- 
termined using a Feldman-Cousins approach 18] . We 
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FIG. 4: Limits on N° L pair production as a function of its 
mass. The limit is for a lifetime of 4 x 10 _11 s. 



improve on the NuTeV limit by several orders of magni- 
tude at long lifetimes and add coverage at lower lifetimes. 
Our limit excludes the interpretation of the NuTeV ex- 
cess as arising from any model with similar N^ L produc- 
tion cross sections and kinematics. 

To summarize, we have presented an analysis sensitive 
to neutral, long-lived particles decaying to pp + X using 
a new technique that expands the capabilities of the DO 
experiment. The background is estimated to be 0.75 ± 
1.1 (stat) ±1.1 (syst) events. The signal region contains 
events and a limit is set. The 95% CL limit for a mass 
of 10 GeV and lifetime of 4 x 10" 11 s is 0.14 pb. This 
result excludes an interpretation of the NuTeV excess of 
dimuon events in a large class of models. 
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